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Abstract. - An obvious determination of the acceleration of gravity g can be deduced from 
the measurement of the velocity of falling atoms using a 7r — 7r pulses sequence of stimulated 
Raman transitions. By using a vertical standing wave to hold atoms against gravity, we expect 
to improve the relative accuracy by increasing the upholding time in the gravity field and to 
minimize the systematic errors induced by inhomogeneous fields, owing to the very small spatial 
amplitude of the atomic center-of-mass wavepacket periodic motion. We also propose to use 
such an experimental setup nearby a Watt balance. By exploiting the g/h (h is the Planck 
constant) dependence of the Bloch frequency, this should provide a way to link a macroscopic 
mass to an atomic mass. 



Introduction. - The dynamics of an atomic wave packet in a periodic potential under the 
influence of a static force has been extensively analyzed using different physical approaches : 
in terms of Wannier-Stark resonance states [1] , Bloch oscillations [2] or macroscopic quantum 
interferences induced by tunnelling due to the external acceleration [3] (for review see [4]). 
An interesting configuration occurs when the external force is induced by the acceleration of 
gravity. In this case the Bloch frequency is equal to vb = ^mr- and depends only on the local 
acceleration of the gravity g, the wavelength of the light A and some fundamental constants. 
This frequency is typically in the range 100 Hz — 2000 Hz, and its measurement allows the 
determination of g. Previous experiments have already been realized using the dynamics of 
BEC [3,5] or degenerated Fermi gas [5] in vertical 1-D optical lattice. Kasevich's group has 
observed the interference between macroscopic quantum states of BEC atoms confined in a 
vertical array of optical traps. This interference arises from the tunneling induced by the 
acceleration of gravity and appears as a train of falling atomic pulses. The acceleration of the 
gravity g was determined by measuring the spatial period of the pulses train. In Ingusio's 
group [5] the Bloch period is straightforwardly deduced from the evolution of the momentum 
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in the trap by adiabatically releasing the cloud from the lattice. In both experiments the 
detection is performed by imaging the falling atomic cloud using absorption imaging techniques 
and the uncertainty on the g determination did not exceed 10 -4 dominated by the imaging 
system. 

Our experimental approach is based on the precise determination of the velocity distribu- 
tion of atoms along the vertical axis using Doppler-sensitive Raman transitions. An obvious 
determination of the acceleration of the gravity is possible by measuring the atomic velocity 
variation after a given falling time T. This is performed by applying a 7r — 7r pulses sequence 
with a spacing time T: the first pulse defines an initial velocity by selecting a narrow velocity 
class from an ultracold atomic sample. Atoms are then in a well defined internal state. The 
second pulse measures the final velocity distribution of the atoms after the fall by transferring 
a resonant velocity slice to another internal state. This so-called velocity sensor allows us 
to locate the center of the velocity distribution with high accuracy and is now limited by 
the experimental setup platform's vibrations [6] . We could substantially improve the relative 
uncertainty on the measurement of g by increasing the falling time T, but this parameter is 
swiftly limited by the dimension of the vacuum chamber. In our experiment we suggest to 
hold the atoms against gravity by applying between the two Raman pulses a far resonant 
standing wave, during an interrogation time TbiocH (fig- Ha). In the non dissipative case, 
atoms fall until they absorb a photon from the upward wave and emit a stimulated photon 
in the downward wave. The atoms make a succession of A transitions (fig. Hb) inducing a 
momentum exchange of 2hk [k is a wave vector of the standing wave) with the cycling fre- 
quency vb — ^gr- • This evolution is equivalent to the dynamics of the Bloch oscillations [7] . 
In a previous work [8] we have measured a transfer efficiency of 99.5% per cycle. This result 
promises a large number of Bloch oscillations. Another particular interest of our method is the 
small spatial amplitude of the atomic wavepacket motion, during the oscillations ( 1 ), therefore 
the selection (first Raman pulse) and measurement (second Raman pulse) are done in a small 
volume allowing a better control of systematic effects arising from inhomogeneous fields. In 
this letter we investigate the possibility to make a high precise measurement of the acceler- 
ation of gravity g using Bloch oscillations of cold atoms in a vertical standing wave. Such 
accurate measurements have important repercussions on geophysical applications including 
earthquake predictions, locating oil and studies of the global warming. We first describe our 
measurement method, then we present a preliminary measurement of g, and finally we discuss 
the signal losses observed when we increase the number of Bloch oscillations. 

Experimental Set-up. - The main experimental apparatus has already been described 
in reference [8]. Briefly, 87 Rb atoms are captured, from a background vapor, in a <j + — a~ 
configuration magneto-optical trap (MOT). The trapping magnetic field is switched off and 
the atoms are cooled to about 3 fiK in an optical molasses. After the cooling process, we 
apply a bias field of ~ 100 mG. The atoms are then optically pumped into the F — 2, rap = 
ground state. The determination of the velocity distribution is performed using a tt — it 
pulses sequence of two vertical counter-propagating laser beams (Raman beams): the first 
pulse with a fixed frequency v se i, transfers atoms from 55*1/2, \F — 2, mp = 0) state to 55*1/2, 
\F = l,mp = 0) state, into a velocity class of about « r /15 centered around {\v se i ect /2) — v r 
where A is the laser wavelength and v r is the recoil velocity. To push away the atoms remaining 
in the ground state F=2, we apply after the first 7r-pulse, a laser beam resonant with the 
55*1/2 (F = 2) to 5Pi/ 2 (F — 3) cycling transition. Atoms in the state F=l fall under the 

( 1 )The amplitude of the oscillation of atomic wavepacket is given by Az = A n /2\F\, where A n is the energy 
width of the nth band [2]. In our experiment only the fundamental band is considered Ao < E r , F = mg 
Then Az < 1/im which is very small compared to the size of the atomic cloud (1 mm). 



P. CLADE et al: MEASUREMENT OF THE LOCAL ACCELERATION OF GRAVITY USING BLOCH OSCILLATIONS3 



a) 

71 - Pulse n - Pulse 
Vertical standing wave 


b) A 

uVW / \ WW 




< ► 

Tbioch * 


* ► 

T 


-/ik fik *p 



Fig. 1 - a). Experimental pulses sequence, b) A narrow velocity class is selected by the first Raman 
pulse. When the vertical standing wave is switched on, atoms fall until they are resonant with the 
A transition: they absorb a photon from the upward wave and emit a stimulated photon in the 
downward wave. This induces a momentum exchange of 2Tik. 



acceleration of gravity during T. We then perform the final velocity measurement using the 
second Raman 7r-pulse, whose frequency is v meas . The population transfer from the hyperfine 
state F — 1 to the hyperfine state F — 2 due to the second Raman pulse is maximal when 
2n(i^ se i—i' meas ) — gxTx || (ki— k 2 ) ||, where ki, k 2 are the wave vectors of the Raman beams. 
The populations (F = 1 and F — 2) are measured separately by using the one-dimensional 
time of flight technique developed for atomic clocks and depicted in [9]. To plot the final 
velocity distribution we repeat this procedure by scanning the Raman beam frequency v meas 
of the second pulse. 

The two Raman beams are generated using two diode lasers injected by two extended- 
cavity diode lasers (ECLs). To drive the velocity-sensitive Raman transition, the frequency 
difference of the master lasers must be precisely resonant with 87 Rb ground-state hyperfine 
transition (~ 6.8 GHz). The frequency of one ECL is stabilized on a high stable Zerodur 
Fabry-Perot (ZFP) cavity. The second ECL is then phase- locked to the other using the beat 
note technique (see fig. 0b). The very stable RF source is performed by mixing the 62th 
harmonic of a 100 MHZ quartz oscillator with different digital synthesizers (SRS DS345). 
They are used to tune finely the frequency of the RF source. A YIG oscillator is phase-locked 
on to the central line of the source in order to reject completely the residual sidebands of the 
different mixings. A multiplexer switches between two synthesizers to generate the frequency 
offset for the velocity selection or measurement. The frequency of a third synthesizer is 
linearly swept during the selection and the measurement pulses to compensate the Doppler 
shift during the fall of the atoms (fig.EJa). The time interval between the two Raman pulses 
is precisely defined by the delay between the triggering signal of the two frequency sweeps. 
Each laser beam passes through an acousto-optic modulator (~ 80 MHz) for timing (switch 
on and off) and intensity control. The two beams have linear orthogonal polarizations and 
are coupled into the same polarization maintaining optical fiber. The pair of Raman beams 
is sent through the vacuum cell. The counter propagating configuration is achieved using a 
polarizing beam-splitter cube and an horizontal retroreflection mirror placed above the exit 
window of the cell. The standing wave used to create the 1-D optical lattice is generated by 
a ThSapphire laser, whose frequency is stabilized on the same highly stable ZFP cavity. This 
laser beam is splitted in two parts. To perform the timing sequence, each one passes through 
an acousto-optic modulator to control its intensity and frequency. The beams are detuned by 
260 GHz from the 55i/2 — 5P3/2 resonance line to avoid spontaneous emission. With these 

laser parameters, the optical potential depth Uq equals to 2.7 En (En = is the recoil 
energy). For this value, when the external acceleration is due to gravity, the transfer of atoms 
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Fig. 2 - a). Temporal variation of the frequency difference of the two Raman beams, b). The synthel 
and the synthe2 allow us to switch between the selection and a measurement steps. To compensate 
the Doppler shift during the falling of the atoms the frequency of the synthes3 is swept linearly during 
the Raman pulses. 

Fig. 3 - The final distribution velocity of atoms after a free fall of 30 ms. The center of the velocity 
distribution is located with an uncertainty of 3.2 Hz, corresponding to a relative uncertainty of 
4 x 10~ 6 on the measurement of the acceleration of the gravity . 



to the higher bands remains insignificant for several periods of Bloch oscillations. 

Results. - In a first experiment, we determine g by measuring the atomic velocity varia- 
tion after the free fall of atoms during 30 ms, using the tt — tt Raman pulses sequence described 
previously. The typical final velocity distribution is shown in (fig. 0J). The center of this dis- 
tribution is located with an uncertainty of 3.2 Hz (corresponding to v r /5000) in an average 
time of 20 ran, allowing a measurement of g with a relative uncertainty of 4 x I0~ 6 . This 
uncertainty is limited by many systematic errors. These errors may occur due to the vibration 
noise of the retroreflecting mirror [6] , the fluctuation of the number of detected atoms and the 
atomic motion between the two pulses (effect of inhomogeneous fields). 

In a second experiment, we apply between the two 7r-pulses a standing wave during an 
interrogation time Tsioch- We then study the evolution of the final momentum distribution 
by changing Tsioch- Before analyzing the experimental results, we briefly recall the relevant 
results of Bloch's theory. The energy spectrum of the particle presents a band structure 
(indexed by n) arising from the periodicity of the potential (optical lattice with period d = 
A/2). The corresponding eigenenergies E n (q) and the eigenstates \n,q) (Bloch states) are 
periodic functions of the continuous quasi-momentum q, with a period 2k = 2n/d. The quasi- 
momentum q is conventionally restricted to the first Brillouin zone ] — ir/d,ir/d]. If we apply 
a constant force F, sufficiently weak in order to avoid interband transitions, a given Bloch 
state \n, q(0)} evolves (up to a phase factor) into the state \n,q(t)) according to 

q(t) = q{0) + 2k— (mod 2n/d) (I) 
tb 

When the atoms are only submitted to the gravity force, the Bloch period tb is given 
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by This period corresponds to the time required for the quasi-momentum to scan a 

full Brillouin zone. In our experiment, first we prepare Bloch states around q = (in lattice 
frame) at the bottom of the fundamental energy band (n = 0) by turning on adiabatically the 
standing wave (rise time of 300 /is): this avoids a transfer of population into the higher energy 
band. We point out that just before turning on the Bloch potential, selected atoms reach a 
mean velocity of about 10 v r . In order to compensate this velocity drift, the upward beam's 
frequency is shifted by ~ 150 kHz (in the laboratory frame the standing wave is then moving 
with a constant velocity of about 10 v r ). To use a pure standing wave, we should launch 
atoms in ballistic atomic-fountain trajectories either from a moving molasses [10] or with 
Bloch oscillations [8], and turn on the Bloch potential when they reach their summit. After 
time Tb we suddenly switch off the optical potential and we measure the final momentum 
distribution in the first Brillouin zone. In (fig. 01a) we report the center of the measured peak 
as a function of the holding time Tg . The observed sawtooth shape is the signature of Bloch 
oscillations (eq.[IJ. We observe more than 60 Bloch periods corresponding to the longest lived 
Bloch oscillator observed in bosonic systems. To determine the Bloch period tb, we measure 
the time interval between the centers of the two extreme slopes of sawtooth. We extract the 
value of tb by dividing this time interval by the number of periods. This measurement leads 
to a determination of the local acceleration of gravity with a relative uncertainty of 1.1 10 -6 
as jjj ratio [11] and the wavelength A are known with a better accuracy. The linear fit of 
the experimental data in (fig. |Ua) is performed by fixing the value of the recoil velocity. In 
fig- Gib we present the residuals of the fit; they increase when we move away from the center 
of the Brillouin zone. That reveals that there is a difference between the quasi-momentum 
and the measured momentum of about 10~ 3 x hk. The momentum spectrum ^(^(p) can be 
calculated by projecting the Bloch state |0, q) onto plane-wave components \p) (measurement 
basis). Using analytic properties of the Bloch wavefunctions one obtains 

*o,,(p) = J^*o(p)x$>(p-<Z-27rZ/d), leZ (2) 

I 

where $>o(p) is the Fourier transform of the Wannier function of the fundamental energy 
band [12]: the momentum spectrum is composed of peaks separated by the reciprocal-lattice 
vector, 2k — 2w/d, with an amplitude given by the envelope function $o(p)- This function is 
not constant along the width of the selected velocity class. Thus the center of the measured 
momentum distribution is shifted with respect to the quasi-momentum q(t) (eqQ by a factor 
depending on the momentum spread. The data analysis allowing the determination of the 
Bloch period is performed in order to reduce substantially the systematic error induced by 
this effect. 

When we increase the interrogation time of the standing wave up to 100 ms, the signal is 
significantly degraded (loss rate becomes larger than 50%). In order to understand the origin 
of this losses, we have measured, for a given laser intensity, the losses per oscillation versus 
the detuning A, relative toF = 2^F = 3 transition (fig. EJ). These losses are obtained 
by comparing the number of atoms measured after N oscillations to those measured after 
N = 10 oscillations. The choice of the parameter A allows us to estimate the effects of the 
spontaneous photon scattering and the interband transitions which depend both on A (for 
interband transition see [7]). We observe that the losses exceed 3%; they are more important 
than the 0.5% rate losses measured using an accelerated standing wave [8]. This results from 
the slowness of the Bloch oscillations, since in a vertical standing wave the Bloch period is 
^1.2 ms, when in reference [8] the period was only ~ 0.1 ms. 

For the experimental values of the standing wave parameters, the losses induced by the 
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Fig. 4 - a. The center of the final velocity distribution versus the duration of the standing wave. The 
dot represent the experimental data and the line the least-square fit performed by fixing the recoil 
velocity, b) The residuals of the fit. 

Fig. 5 - The losses per oscillation versus detuning A. These losses are obtained by comparing the 
number of atoms measured after TV oscillations to those measured after N = 10 oscillations. The 
relevant advantage of this presentation is to take into account only the losses during the Bloch 
oscillations process. Dot: the experimental data. The theory: (Dotted line). Solid line: theory 
including the interparticle collisions. 



spontaneous emission and the interband transitions in the weak binding limit do not match 
with the experimental data (dotted line in (fig. 0). Performing a least-square fit based on 
this model and including collisions with the residual Rb vapor, we extract a characteristic 
time constant of the damping due to the collisions of about 70 ms. This value corresponds 
to the lifetime of the molasses for the residual vapor pressure in the cell. Therefore we think, 
that the number of Bloch oscillations, in our experiment could be increased by reducing the 
pressure in the vacuum chamber. 

Conclusion and prospects. - We have described an experimental method to measure the 
vertical velocity distribution of atoms by using a n — it pulses sequence of Doppler-sensitive 
Raman transitions. We have performed a preliminary determination of the local acceleration 
of gravity with a relative accuracy of 10 -6 by measuring the Bloch period. We have also 
demonstrated that the number of Bloch oscillations is not yet limited by either the interband 
transitions or the spontaneous emission, but only by collisions with the background atomic 
vapor. To overcome this limit, we are now building a new ultrahigh vacuum chamber where 
the magneto-optical trap will be loaded by an atomic slow beam. In order to improve the 
accuracy of the velocity measurement a vibration-isolation system is also in implementation. 
These improvements should allow us to take a better benefit of the Bloch oscillations. 

An attractive possibility consists to replace the 7r — n velocity measurement by a two 
7r/2 Ramsey-Borde sequence. We obtain then a ir/2 — 7r/2-Bloch oscillation- tt/2 — ir/2 atom- 
interferometer. Comparing this scheme to the 7r/2 — 7r — 7r/2 atom-interferometer used in 
gravimetry [13, 14], where the pulses spacing time is limited by the effects related to the 
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spatial position, the advantage would be to increase the measurement time thanks to Bloch 
oscillations for a similar phase difference between the two paths. This way, we could reduce 
significantly the uncertainty of the gravity interferometric measurement. Finally, we suggest 
to use such experiment nearby a Watt balance site [15, 16]. In the dynamics mode of the 
balance the relation which equates the mechanical power and the electrical power is given 
by [16]: 

Mgv = CFjF'jh (3) 

where M is the standard mass, v is the velocity of the vertical moving coil, h the Planck 
constant, C represents a dimensionless constant and Fj, F'j denote the frequencies applied to 
a Josephson device. Using the g/h dependence of the Bloch period we obtain 

M CFjF'j A 1 

TO 2 Vq V 

In conclusion, associating our Bloch oscillations experiment to a Watt Balance could be used 
to link a macroscopic mass to an atomic mass. 

* * * 
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